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Summary 


A theoretical method is established for determining the aero- 
dynamic ch;iracteristics of over-wing-blowing configurations. The 
method accounts for both jet entrainment and jet interaction effects 
because of the differences in freestream and jet dynamic pressures 
and Mach numbers. The predicted lift increments agree well with 
available data. It is shown that the lift is under predicted with 
entrainment effect alone when the jet is close to the wing surface. 



1. Introduction 


The over— wing engine installation has recently i been investigated 
as an alterative of installing high by-pass ratio turbofan engines. 

With the conventional under— wing installation.} these big engines may 
have ground clearance problem and adverse aerodynamic effects on the wing. 

In addition to solving the ground clearance problem, the over-T^ing_ in*^ 
stallation has the advantages of reducing jet noise because of Mng shield- 
ing effect and improving the aerodynamic performance. In some experimen- 
tal investigations} Bagley and his associates tested a model with rectangu- 
lar jet nozzle at high subsonic freestream Mach numbers [Refs. 1 and 2] . 
Similarly, extensive experimental study of various configurations has been 
made at VFW-Fokker in connection with the VFW-614 project [Ref. 3]. In 
the test configurations mentioned above, the nacelles were installed d.Lrect- 
ly above the wing. Therefore, the nacelle-pylon solid-body interference on 
the wing flow will be significant. On the other hand, Seidel [Ref. 4] and 
Putnam [Ref. 5] measured wing aerodynamic characteristics with the jet exit 
ahead of the wing leading edge. Seidel also allowed the blowing jet to be 
inclined relative to the freestream. From these limited experiments, the 
favorable aerodynamic effects of the over-wing blowing jets have been con- 
firmed. 

In general, when a turbulent jet is blowing above the wing, the wing 
upwash will be increased due to jet entrainment. The jet entrainment ef- 
fect on the wing can be simulated by using sink distribution albng the jet 
axis. For this purpose. Squire and Trouncer [Ref. 6] derived a theoretical 
sink strength distribution for incompressible, non-heated jets. Thxs theory 
has been used to .estimate the jet effects on the horizontal taxi [Ref. 7] . 
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Ribner extended Squire and Trouncer's analysis to heated but incompressible 
jets [Ref. 8], In these early investigations;,, , the primary .ob^^^ , 

been the study of jet entrainment eff^tts pn the tailv ,. ph t^^^ 

Putnam [Ref. 9] has recently applied Squire and Trouncer^s theory to the 
investigation of the wing aerodynamic characteristics due to over-wing 
blowing jets at moderate to high subsonic Mach nimhersl; Because the theory 
is applicable only to incompressible flow, Putnam used equivalent incom- 
pressible jet parameters in his applications. Line sink distribution was 
used in conjunction with the vortex-lattice method. Another approach used 
at VT'ff— Fokker was to use sink panels over the jet surface [Ref , 3] . Both 
the sink strength and the jet bomdary for the velocity ratio V /V. of 1/7 
were based on the static jet properties, so that the entrainment used in ■ 
the computation is expected to be higher than the actual magnitude. In 
both Investigations, it has been indicated that when the jet is close to 
the wing surface, the methods xrould be inadequate. In fact, the predicted 
lift due to entrainment alone would be too loxif as compared mth the experi- 
ments [Ref. 9]. Therefore, additional lift mechanism must be examined. 

In the present investigation, the lift induced by the over-wing 
blowing jets will be shoxm. to be attributable to the jet entrainment 


and the inviscid jet-wing interaction due to the higher dynamic pressure 
in the jet. The latter effect is particularly important x^hen the jet is 
close to the wing surface. To extend the capability of existing methods 
for computing jet entrainment to compressible heated jets, Kleinstein's 
theory for static compressible jets [Ref. 10] is extended to the case 
of external stream, with the simplicity of his theory being retained. 

The wing-jet interaction theory xcith Mach number nonuniformity reported 
earlier [Ref. 11] will then be used to treat the interaction betx?een the 
x^ing and a rectangular or a round jet. 
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2. List of Symbols 
AR wing aspect ratio 

bi/i nondimensional turbulent jet baif width defined in Eq. A. 21 
c local chord length, m (ft) 

c reference chord length, ro (ft) 

C leading-edge singularity pat imeter (see Eq. (27)) 

sectional induced drag coefficient 
“i 

Cj flap-chord length, m (ft) 

sectional lift coefficient 

AC- difference in lift coefficients with jet on and off 

Lj 

c sectional pitching-moment coefficient about y~axis 

m 

pressure coefficient, or specific heat at constant pressure 
c^ sectional leading-edge thrust coefficient 

d^ initial jet diameter, m (ft) 

^ ^ 
e unit vector tangent to jet path and is taken to be i in the 

shallow jet approximation 
E(x) entrainment function defined in Eq. A. 65 
F]^ a function defined by Eq. A. 11 

F^ a function defined by Eq. A. 12 

F 3 a function defined by Eq. A.38 

H stagnation enthalpy 

i,j,k unit vector along x-, y- and z-axes, respectively 

kj a constant in Ferri's turbulent eddy viscosity model 
(see Eq, A. 20) 



ni mass flow rate,Hg/sec (slugs/sec) 

m’ slope 

M Mach, number, or momentum 

N number of chordwise integration points 

[N] . normal velocity influence-coefficient matrix 

n unit vector normal to the jet surface 

n,s jet axis system with n normal to the surface and 
tangential to it 

p free stream static pressure, N/ra^ (Ib/ft^) 

CO 

p . jet stagnation pressure, N/m^ (Ib/ft^) 
t j j 

P^(x) a function defined by Eq. A. 16 
PjCx) a function defined by Eq. A. 10 


? 2 (x) a function defined by Eq. A. 13 


q entrained volume flow rate, m^/sec (ft^/sec) 

r radial coordinate, m (ft) 

r^ equivalent jet radius, m (ft) 

initial jet radius, m(ft) 


R 

R. 

J 

S 

[S] 

t 

T 

T 


average gas constant for the air-jet mixture 
gas constant for the jet 

wing area, m^ (ft^) 

tangential velocity influence-coefficient matrix 
airfoil thickness, m (ft) 
temperature, (°R) 


jet stagnation temperature, °K (“R) 



free stream stagnation temperature, K ( R) 


velocity in x-direction, m/sec (ft/sec) 

= u/u^ 

equivalent jet velocity, ra/sec (ft/sec) 
velocity, m/sec (ft/sec) 

unperturbed velocity vector, m/sec (ft/sec) 

perturbed velocity vector, m/sec (ft/sec) 

jet— entrained— flow velocity vector, m/sec (ft/sec) 

= x/r 

o 

nondimensional turbulent jet potential core length based on r^ 


jet exit X- and z-coordinates 

wing— fixed rectangular coordinate with positive x-axis along 
axis of symmetry pointing downstream, positive y-axis pointing 
to right and positive z-axis pointing upward, m (ft) 

coordinate of camber surface, = z^(x,y) , m (ft) 


a angle of attack, deg 

g - (1-M^)^ 

Y nondimensional vortex density, or ratio of specific heats 

6 airfoil trailing-edge half angle, deg 

e 

8 . jet deflection angle, deg 

2 

Cl - r/r^ 

0 = T . /T 

3 e 

e(x) angle of camber slope 

A sweep angle, deg 
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V == Vjv., or n^/n. 

y» =3 V .e/V. .e 

o j 

p density, kg/m^ (slugs/ft^) 

■p. equivalent jet density, kg/m^ (slugs/ft^) 

1 

^ Tiondimensional turbulent transfer coefficient 

^ dimensional perturbation velocity potential, ni? /sec (ft^/, sec) 

^ nondimensional perturbation velocity potential 

nondimensional additional perturbation velocity potential, : 
or stream function 


Subscripts: 
a additional 

c jet centerline 

e external or outer flow 


j 

jj 

£ 

o 

oj 

t 

w 

wa 

Mj 

wo 


jet flow 

jet vortices for jet flow 

leading edge 

outer or external flow 

jet vortices for outer flow 

trailing edge 

wing 

additional wing vortices 
perturbation due to wing in jet flow 
perturbation due to wing in outer flow 
free stream 
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3. Mathenjatical S’ornmlation =; - 

3.1 Inviscid Interaction 

In general, a jet mil expand after leaving the nozzle and its 
axial velocity will decay dowstream. The velocity distribution over 
any jet section is nonuniform, with maximum value at the center, for 
simplicity, an equivalent jet section mth uniform velocity distri- 
bution, to be defined later, will be used in the ptesent formulation. 
Since the jet flow has higher dynamic pressure than the freestream, 
any disturbances created outside the jet mil be reflected and dif- 
fracted at the jet surface [Ref. 12]. The disturbances in the present 
applications are mainly due to the presence of the wing in the free— 
stream and the jet-entrained flow. This inviscid wing-jet interaction 


is the subject of fomulation in the present section. 

The perturbations of the flow field mentioned above are determined 
by satisfying the boundary conditions. On the jet surface, it is re- 
quired that the slopes, of the streamlines on both sides of the surface 
are the same. When linearized, this condition can be e^qpfessed as 


[Ref. 11] . 

Ik . . w 

One additional jet surface boundary condition is the static pressure 
continuity, which in linearized version can be written as 

In the above, e is a unit vector tsaigeut to the jet path and is taken 

to be I in the shalloW jet approximation, and ri is the mit vector - 
normal to the surface. VV is the: unperturbed jet velocity and is ■ 
the external velocity vector before the ^H.ng is intfoduced into the. 
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flow field. 


With fhe wing being assumed to be entirely in the outer-flow 
region, the wing tangency condition has the conventional form: 

= - V-S <«. 

For numerical convenience, Eqs. (l)-(3), xd.ll first be rex^ritten 
in. a npndimensional form. For this p.urpose, define the nondimensional 
perturbation potentials ' arid such that 

Introducing Eqs. (4)-(5) into Eqs. (1)- (3) , it is Obtained xdth so^^ 
simplification that ; ^ 





(6) 

>1 

■vn 


-as 

- 1 

Li ')*- zi. =; 0 

(7) 

it 

V • 

¥|0 


(8) 

■ z 

“ V ^ 

e S 

T 
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(9) 


X7here 


Let V”. be the jet-entrained flox? and the uhifoi^ frOe&trea% 
vector. Then 

V - V +■ V- (10) 

It is assumed that the jet-entrained flox? does not affect the jet stream 
surface condition (Eq.‘ 6) xn the ■ interactxon rormulation. ■ It . xhf luences 
only the x/ing tangency condition through Eq. (8). Since the problem is 
iiriear, it is convenient to decompose it into the xdng-alone case and 
the additional jet— induced perturbation. Therefore, let 



whera are the wing-alone perturbation poten- 

tials at respective Mach numbers. Using Eqs. (11)-(12) , Eqs. (6)-(8) 


become 

'T>n 




•T>yy 


3- 


_ .j. 


'4'^ 




-o-n 


■3 S 
•3 z; 


t(/*'7 ^ 






3S 


'Vro 
3S 


= o 


(13) 

(14) 

(15) 


where and satisfy the conventional tangency condition with 

jet entrainment: 

- +.»=< - 

<z<z \U,«soC 

- +a„.i - >• * - (17) 

'iZ 

Because of Mach number nonuniiormity , two vortex sheets 'ire 
introduced on the jet surface, one with strength to account for 

the outer fl ow perturbation and the other with strength to ac- 
count for the jet flow perturbation, let be the additional wing 

vortex strength. Then Eqs. (13)— (15). can be reduced to ,[R^« 111 

5 _ g.- j c t CIS) 


V^‘ e 


f (.“O'c - T0“'r[ 1 


7-1^): + 




fr; 


09) 


where [K] and [S] arc the influence coefficient matrices for induced 
normal and tangential velpcity. cqmpon^ts , with the , fir^^ 
denoting the control point locations and the seqond f or the influencing 
vortex locations. When the jet is close enough to the wing, it may 
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attach to the \d.ng upper surface due to the Coanda effect. As a 
result, the jet may deflect at the wing trailing edge. This jet 
flap effect can be included in accordance with method described in 
Ref. 11. Both [N] and [S] matrices are formulated with the Quasi 
Vortex— Lattice Method (Quasi VLM) (Ref. 13). Since the Quasi VLM 
is one type of the discrete— element methods, care must be taken in 
evaluating [S] matrices for better accuracy (Ref. 11). 

Once Eqs. (18)- (20) have been solved, the total td.ng vortex 
strength is given by 

/ = y d- / ( 21 ) 

The sectional wing aerodynamic characteristics can be computed by 

properly resolving the pressure force acting on the camber surface 

into appropriate directions. It follows that 

I 

3 V.-?) tu:) 


^ ~ ^ tos(^o4- Bij} (23) 




a «jso<. r * 


c c 


C - 

3 4.(^5 


where 


SC»)= Wc tt ) 


(24) 


(25) 


The sectional leading-edge thrust coefficient c^ is computed 
as (Ref. 13) 


7r 




C !- Mo tos^Ac) ~ 


* X ujs Ajj 

where the leading-edge singularity parameter C is defined as 
C! = litr, si-n© 

o 


(26) 


(27) 
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According to the Quasi VLM, C can be determined by summing the total 
induced normal velocity at the leading edge and subtracting the right- 
hand side of Eq. (17) or (20), depending on whether the effects of 
wing alone or the jet interaction is being considered. If the induced 
normal velocity at the ith leading edge control point is denoted by a 
suVscript i, then the singularity parameter Cj in the wing alone case 
with flap deflection can be computed from the expression; 

C, \ - 5 ^ -tano< (28) 

Similarly, the singularity parameter, due to the jet-induced per- 
turbation can be computed from Eq. (29): 




In the above, N, is the chord^-Tise number of vortex elements on 
the airfoil section, The integrals in Eqs. (22)-(24) are computed with 
the midpoint trapezoidal rule after the 0 -transformation. The over- 
all aerodynamic characteristics of the wing are then determined by 
spanwise integration of the sectional characteristics (Ref. 11). 

3.2 The Jet Entrainment Model 

The jet entrainment affects the wing aerodynamic characteristics 
by creating additional upwash on the wing. This additional upwash is 
rejsresented by the last term in Eq, (16) or (17). The ability of pre- 
dicting the jet-entrained flow in subsonic compressible flow is im- 
portant in the present investigation. Although it is possible to 
solve the jet problem by numerically D.itegrating the governing equa- 
tions (Ref. 14), such an approach would greatly increase the total com- 
puting time. With the idea of having a method being simple, general 
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and accurate in mind, it appears that Kleinstein's theory (Ref. 10) 
is the best choice. However, Kleins tein presented his theory for 
free jets only. In Appendix A, a method is described to extend 
Kleinstein’s theory to compressible jets in an external stream, with 
the primary end results being the prediction of jet entrainment. 

As shown in Appendix A, the entrained mass flow rate is computed 
as (See Eq. A. 42). Let the entrained volume flow rate be de- 

noted by q: 




I -aTn 


(30) 


The effects of the jet entrainment on the nearby lifting surfaces can 
be simulated by a line sink distribution on the jet axis mth strength 
equal to the entrained volume flow rate. The induced velocity com- 
ponents due to this sink distribution can be computed as (Ref. 15) 
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(31) 




■ax 


(32) 


Instead of using q directly, it is convenient to use the nondimeu 


sicnal entrainment function E(x) defined as (Ref. 16) 
dr ' 
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(33) 


Hence 


(34) 


In applications, Eqs. (31) and (32) can be integrated by dividing 
the jet axis into segments within each of which E(^) is assumed to 
vary linearly. Let 

ECi)=^ a, + (35) 
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in the ith interval. Then 


J. 

•ff 




J_ :± — %r ■?- V 

Similarly, 


-i- ^ 


y^ r (ft^-i-ki)c^-s)^j 


w^ro 


The integrals in Eqs. (36) and (37) can be exactly integrated to give 


where 


(a^-f biS3 <^5 a^+bix^ , ^ 

i cu--!)Vf'y-'- " ■ 

<;».<■ Uj)Cx-I? Aj = (a, + b;x)G..,-l>;S,; 




X - X; 


7C - 




Cx- 


Note that only -component will contribute to an upwash on the xjing. 
To find -ft in Eqs. (16) and (17), this -component, which is 

always pointed toward the jet axis, must be decomposed at the \'ri.ng 
plane to produce the upwash- The x-induced velocity is small in 
general and is neglected in the present applications. 

3.3 An Equivalent Inviscid Jet in the Interaction Theory 

As described in Appendix A, the velocity and temperature distri- 


tx- 


butions in the turbulent jet vary not only in the axial direction but 


14 



also radially. The propagation of any disturbances in such a nonhomo- 
geneous region would be governed by a partial differential equation 

non-constant coefficients. To find a solution of such an equa- 
tion would be difficult. To make the solution tractable, an equiva- 
lent homogeneous jet must be defined. In the present analysis, an 
equivalent jet with constant cross section and jet properties is as- 
sumed to be one which satisifes the conservation of mass, momentum 
and heat content and is evaluated at the midchord of the Tring if 
the jet exit is ahead of the leading edge. The use of an equivalent 
jet in the wing-slipstream interaction problem has been shown to yield 
reasonably good results (Ref. 16). The detail of determining the 
equivalent jet properties is described in Appendix B. 

Let be the equivalent jet radius. It is assumed that if the 
turbulent jet as defined in Appendix A intersects the wing surface and 
follows it to deflect at the trailing edge (to be discussed in more 
detail later), an upper-surf ace-bloTTing rectangular jet with lateral 
extent equal to a r. can be used for interaction computation. The 
main reason behind this assumption is that a jet washing the wing sur- 
face most likely has a cross section of half ellipse rather than a circle 
(Ref. 17). A rectangular jet would be most convenient to use. In this 
case, the situation would be that of upper-surface blowing, so that the 
theory of Ref. 11 can be directly used. 

On the other hand, if the jet does not wash the mng, or if the 
jet is strong enough to go straight, the circular shape for the equiva- 
lent jet is retained. In the computation, the circular jet is in turn 
approximated by an inscribed polygon. This is shown in Fig. 1. To see 
hov 7 wall this polygon can approximate the circle, a 2-D case is con- 
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structed. In the 2-D case, the exact solution can be obtained by 
Kohing's image method(Ref. 18). According to the image method, a 
vortex of ilnit strength situated outside of the circular jet will in- 
duce, due to interaction, additional perturbation in the outside flov 
field which can be represented by an image vortex A of strength -“Cs*'-0/Cs'+i) 
at the inverse point and a vortex B of strength at the center, 

where . Tliis is described in Fig 2 . 

In the present model, two vortex sheets are used on the jet sur- 
face, as described in Section 3.1. Through numerical experimentation, 
it was found that one chord length of jet each ahead of the leading 
edge and behind the trailing edge would be sufficient for numerical 
convergence. Seven vortex elements for each streamwise vortex strip 
of the jet ahead of the leading edge and only one element for the jet 
behind tlie trailing edge were used, in order to reduce the computing 
time. It appeared that the trailing jet had negligible effect on 
the induced downwash near the leading edge. Therefore, in the follow- 
ing comparison, the additional induced downwash near the leading edge 
is plotted and compared with the results by the image method. This 
is done in Figs. 3(a) and 3(b) for two different originating vortex 
locations. As shovm, the present approximation appears to be reasonably 
good. Further improvement of the accuracy seems to be limited by the 
accuracy of sidewash computation which is needed in calculating the 
induced normal velocity on the jet surface. 



Figure 1 GediSetric Approximation 
for a Circular Jet 



16 






CO 


Figure 3t)"Concluded 




4. Numerical Results and Discussions 


4.1 Some Observations of Falk's Experimental Data 

For convenience, the jet-wing geometry parameters are defined in 
Fig. 4. 



Before applying the above theory to any configurations, it is 
necessary to know what the jet deflection angle is when it washes the 
wing. The deflection would occur even if the flap angle is zero. This 
is because the jet has the tendency to follow the airfoil upper surface 
and leaves the trailing edge at an angle equal to or less than half 
the airfoil trailing-edge angle. However, since the jet is thick, its 
deflection mil not be expected to be at full trailing edge angle at 
all times. Intuitively, at a given jet exit location, if the jet veloc- 
ity is high relative to the frees treara, the jet will have the tendency 
to go straight. Therefore, the deflection angle will be small. On the 
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other hand, the deflection angle will be larger if the jet axis is 
moved toward the mng surface. Furthermore, if the jet exit location 
is moved for^'rard away from the leading edge, the measured lift seems 
to be increased according to Falk’s data. This is illustrated in 
Figure 5. Theoretically, the jet entrainment would be lower near 
the wing region when the jet exit is farther away from the leading edge 
(i.e. , 2 Cj/d^ = -6). However, the measured lift increment is higher 
with ^ 1* This seems to imply that the jet will deflect at a 

larger angle. In summary, the jet deflection angle for the thick jet 
under consideration is a function of /a , x^/d^ and Zy Without deter- 
mining its correct variation, it would not be possible to predict the 
Interaction effects. Since a theoretical analysis of this problem is 
beyond the scope of the present investigation, empirical relations would 
be useful alternatives. The effects of w and x./d can be combined by 

JO ^ 

using the equivalent velocity ratio JX (i.e., u /u.) defined in Appendix 

® J 

B. The follovjing procedures are then used to derive such empirical 
relations. 

(1) Three /A ’s (0.333, 0.4 and 0.5405) in Falk’s experiments for 
a rectangular Wing of AR = 2 at “ “1 ^j^*^o ~ ^ used for 

correlation. The deflection angles are taken to be those which would 
produce the correct ^ If is the maximum possible deflection 

angle, then the actual jet deflection angle obtained through the above 
correlation using Jj, computed at the ^ri.ng midchord can be written as 
S’j- = -f 5'e (42) 

where 

+ 35.737l> - S'.l i-f? , ^ (43) 


20 




Figure. 5 Effects of Jet Exit Locations 
on. Fallc’s Experimental Lift 
Inc remen ts, y=0.333, ct=0‘’ 


where /I = 0.6339 is the equivalent velocity ratio for /U = 0.5405. For 
Jx "7 0.6339, linear interpolation between f at ^ = 0.6339 and f = 1 at 
yU. = 1 is used. If f is negative by Eq. (43), f is to be replaced by zero 
and the equivalent jet is assumed to be circular, instead of a rectangular 
one. 

(2) For z,/d * 0.75, full deflection angle and rectangular jets 

1 o 

are assumed irrespective of values. 

(3) Linear interpolation is used for 0.75 < ^ l-O* 

Further assumption is made that even if the jet is deflected, no 

jet reaction forces (i.e. Coanda forces) will be produced. This is based 
on the observation of static data which show that the lift force will 
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still be produced even if the jet is moved away from the wing. There- 
fore, the lift is not produced by jet reaction. Also, the fact that the 
jet tends to attach itself to the wing surface is mainly due to the jet 
entrainment, the effect of which has already been included. The static 
lift force may be explained as produced by the jet suction effect (or 
entrainment effect) which is called "Lagally force" in hydrodynamics 
(p. 215, Ref. 19). The Lagally force is neglected in the following ap-. 

plications. 

4.2 Comparison of Predicted Results x^ith Experiments 

Falk’s experiments represented an early investigation of over-wing- 
blowing aerodynamics. Systematic measurements of lift increments 
due to jet blowing have been made for different jet exit locations varied 
both vertically and horizontally. For his rectangular wing of AR = 2, 
which will be used in the following comparison, the airfoil section was 
a symmetric modified Joukowski airfoil of 16% thickness ratio. By "modi- 
fied," it is interpreted here as a Joukowski airfoil with straightened 
trailing edge by drawing tangents to the surface from the trailing edge, 
instead of the usual cusped shape. The trailing edge half angle can 
therefore be obtained in the following way. According to Glauert (p. 75 
of Ref. 31), the symmetric Joukowski airfoils can be approximately repre- 
sented by the following equations: 

X = — cos 9 (44a) 

2 . 

? “ i: c 1 + &) sin e (44b) 

where c is the chord length. The slope is therefore given by 


The trailing edge is at x = -c/2. The equation for the straight line 
which passes through the trailing edge and is tangent to the airfoil is 
then 


-m' C X + -I;-)- 


(46) 


The tangent point location can be obtained from the following equation. 


from which cos0j= » Hence, m* — —0.0918 from Eq. (45) for t/c 0.16. 


It follows that == 5.25°. 

The predicted results for 3/A.sat two different jet exit x-locations 
are compared with Falk’s experiments in Figure 6. Fr’om this comparison, 
the following points should be noted. Firstly, the present theory xvith 
entrainment alone agrees very well with Putnam’s calculation for/U — 0.333 
and X /d = -1. However, the entrainment -alone theory tends to underesti- 

j o 

mate the lift for z./d 4 1-5, as has been noted by Putnam (Ref. 9). Sec- 

J o 

oudlv for z /d ^1.5, the present results with both interactioi; effects 
JO 


and entrainment agrees reasonably well with experiments. For - 0.25, 

X /d = -1 and z /d =0.75, the prediction appears to be too high, prob- 
j O JO 

ably because it was assumed that the jet was deflected at the full trail— 
ing-edge half angle. Thirdly, as the nacelle is raised upwards, the the- 
oretical method tends to overpredict the lift, even without interaction 
effects. One possible explanation is that the wing might have been af- 
fected by the low dj^naraic— pressure viscous layer from the nacelle which 
is ahead of the leading edge. Such effect has been noted by Stuper in 
his experiments on wing— slipstream interaction [Ref. 32]. 


To test the theory with compressibility effect, Putnam's experiments 
on a swept wing with the engine pressure ratio of 1.9 are compared. The 
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p=0.5405 


Fir,ure 6 Comparison of Predicted Results ^wi tlx Falk, s 
Experiroants . II =0, M.=0 and a=0 







experimental data are sliown on Figure 18 of Ref. 5. Since the velocity 
ratio, temperature ratio and the jet Mach number are needed, they will 
be computed in accordance with isentropic relations. (See Chapter 2 of 





1 + 



X. 


I + 


(50) 

V.. 

= /* = - 

/ 

M,-/ 

s 

(51) 


Using the above relations, it is found that for Pj.^j/p^ “ 1*9, 1.003 

5'»/ ^3 = 0.85914 and = 0.43026. Note that > 1 is allowed in the 

entrainment computation. In the interaction computation, the jet Mach 
number used in the program is given by 



where U: and are equivalent values evaluated at midchord. The equivn 
lent jet Mach number will be much less that one in the present applica- 
tion. The trailing— edge half angle for N.4CA64A006 airfoil u-ed in the 
experiments is found to be 3.6". The results are shotm in Figure 7. 

The agreement with experiments is good. 

Finally, Seidel’s experiments will be examined and compared. In his 
test set-up, the jet exit was at least one chord length ahead of the lead- 
ing-edge. Since the test speed was low, the freestream and jet Mach 
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Figure 7 Comparison of Predicted Results 

with Putnam's Experiments. z^/d^= 

0.75, ./P =1-9 and a=0'’ 

t,i “ 


numbers can be assumed to be zero. The trailing-edge half angle for 
NACA 0010 is found to be 6.67“. The predicted results are compared with 
the measured data in Figure 8, The computed magnitude of A is slightly 
lower than the measurement probably because part of the wing may be im- 
mersed in the jet as indicated by the present turbulent jet theory. 

RTien part of the wing is immersed in the jet, the wing-slipstream inter- 
action theory (Ref, 34) should be used. However, this has not been done 
in the present preliminary investigation. 

5. Cbncl \=;ions 

\ 

A theor itical method has been described for determining the aero- 
dynamic chara -.teristics of over-wing-blowing configurations. The method 
accounts for not only the jet entrainment but also the jet interaction 
effects due to the differences in freestream and jet dynamic pressures 
and Mach numbers. The jet entrainment is determined by a new method 
which is an extension of Kleinstein's free turbulent jet theory. Com- 
parison of the theoretical results with some available data indicates 
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Figure 8 Comparison of Predicted Results with 

Seidel’s Experiments. 

z . / d =1 and a=0 “ 

J o 


reasonably good agrep;. 3 ent. From the examination of some experimental 
data and the theoretical entrainment function, it appears that the en- 
trainment-induced lift can be maximized by placing the wing in the re- 
gion of maximum entrainment. However, any favorable wing- jet interaction 
including the jet flap effect should also be examined for an optimum 
engine location. 
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Appendix A 

An Extension of Kleinstein's Compressible Turbulent Jet Theory 
to the Computation of Jet Entrainment 


A.l. Turbulent Jet Model 


The momentum integral relation for the jet can be written as 
r'^i Y-’- 

J -U..,) (a.i) 


where is the radius of the jet boundary and is the jet radius 
at the exit. For the fully developed region, it may be assumed that 
(Ref. 21, p. 306) 




u = + cw.c-uoci~r"> 


' = Vr. 

Define the nondimensional variables : 


(A. 2) 
(A. 3) 
(A. 4) 


It follows that Eq. (A. 2) can be written as 

Assuming constant pressure mixing, the density can then be written as 

I 


(A. 5) 


(A. 6) 


- _ -RjTv Tj'Ki 


(A.7) 


1?T T<^-R 1+ C^-lXl-i;'-) 

where Rj and R are the gas constants of the exit jet and the jet-air 
mixture, respectively. It follows that Eq. (A.l) can be reduced to 

CA.3, 

where R is the average gas constant for the mixture. 

The integrals Invovled in Eq. (A. 8) can be integrated exactly. How- 

T* 

ever, the results are difficult to evaluate numerically when— ^ - I 

Te 

T 

is small. Since '-I ) is, in general, not large in the fully- 

t e 

developed region, it is convenient to develop the integrand in a 
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power series with respect to ~| ) and integrate the results exactly. 

Te 

In order to accelerate the series convergence. Shanks' Transformation 
(Ref. 22) can be used, l^en the above-nentioned process is applied to 
Eq. (A. 8), it is obtained that (see Appendix C for one example) 


where 






(A. 9) 


(A. 10) 


= d- o.oiGnVi iC+o-opijs-l^^ 

J-f e.8|?| 7 

F =\'_A£tLit2: 

. I + T.ci-5f''.) 

^ 0.06676 4* O-Po'f-S'B.T^ - 0 - Poa-oi-Tl °-Ooo7-S“R 

) + O.gs 7 » 6 71 


(A. 11) 


(A. 12) 


TlCx) = - 1 (A. 13) 

' e 

Eq. (A. 9) gives a relation between the jet radius r^ and the axial 
flow properties. Here, the axial flow properties will be obtained by an 
extension of Kleinstein’s theory (Ref. 10). 

Extension of Kleins tain's theory 

Kleinstein's theory is based on the linearization of the jet flow 
equations on von Mises' plane. The longitudinal velocity ix.»U./uj at 
any x is expressed in terms of the P-function as 

where P* is the normalized P-function and P^ is related to U,t by 


u... + Cl ->w') 

To = I - 


(A. 15) 
(A. 16) 
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(A. 17) 


The 


^ - variable is related to x - ^c/r^ by 




wher® 5^ is the non-dimensional turbulent transfer coefficient. 

The relation between''!' , the stream function, and the radial coordi- 
nate r is given by 

T (A. 18) 




S H 


where 


r= r/n (A. 19) 

According to Ferri's turbulent eddy viscosity model, in tbe 
fully developed region can be assumed to be 

fe. = hy^ C !<. ii-c “ i a) 

where "b,/ is the raidus at which the flow properties satisfy 

To use the above expressions, consider the asymptotic expansions for 
i.e., far downstream. In this region, S’c 

Eq. (A. 21) becomes 


. Hence , 




f =. 




From Eq. (A. 14) 


or. 


Now, 


> + 0-/<) T. = ■L[/x+/- + OV->'P.l 

y> t>*(JL , (A.22) 


tS e T.— i 


C = ± f’’’ '- -i' Ar ^ 4, 




"I* 




S U 






■ir- 




I + 




Cj->^7P 
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But 




by Eq. (A.22). Hence 

T*" =1 f * 

^>'=- te/" 

Uelng Eq. (A.15), Eq. (*-20) cen be written ee. wten | 


(A. 23) 


= c , v'>I A. 

Eq. (A. 24) is va'5lnly Ear dcJcl^T I extrapolate this result » 

the near field. It Is most convenient to use P^Cp. That is, ,5 

Eq. (A.24) is to be replaced by P„. According to Eqs. (A.IO) and (A.15), 

placed by P| C^) • Hence, 

^ ^ n^r St _i±-:feA 

ft - J 7 . ^ L 

For/x=Og Eq. (A. 25) is reduced to Kleinstein’s results for free jet 

fs. - -fe, -/Te , * 

tr a relation between ^ and x can be obtained: 

Integrating Eq. (A. 17), a reiacion u 5 

(A. 26) 


-JE ^ 




,W_JEiSL 

, 4 . 


-V C. 


(A.26a) 


= Cl- e 

For/i=0, Eq. (A.26) gives 

. , Tj -1 - it is obtained that C,- U.J3 

Comparing this equation with Kleinstein , 

which was obtained from experimental correlation. 
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The constant k, will be chosen later for correlation with experimental 
data. Eq. (A. 26) may be regarded as a nonlinear algebraic equation for 
^ (x) . For a given x, the corresponding ^ 5 and therefore U.^ , can be 

obtained by Newton's method. Eq. (A. 26) also shows that the potential 


core length is given by 

C, 


D.SS' 


- 


(A. 27) 


The integral involved in Eq. (A. 26) can be evaluated as follows, 


For large ^ , the integrand becomes j V - . Henc-e ^ 






fnf. 


. . .Su 

Tt 








V i .(A.28) 






where 


i±^ 

1+ 




i-SC* a.^/ 


(A. 29) 


To complete the evaluation of flow properties on the jet axiS; 
it is necessary to determine T - According to Kleins tein s theory. 


the stagnation enthalpy on the axis is given by 

~ WSW - 1 - 


(A. 30) 
(A. 31) 


£ll « 

Tt 

where P is the turbulent Prandtl number. From Eq. (A. 31), it is seen 
r 

that Eq. (A. 26) can be modified to give 




A 



+ 0 . 2 .^ 


(A. 32) 


1 + 


Once H is obtained, T can be obtained from the definition of the 
o c 


stagnation enthalpy: 

H == CpT + -t V." 

It follows that 


(A. 33) 
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He- He _ u - A- ju-e j 

n o 




,Te + ^Uj‘. 

C^Tj + -I: Uj - ( Cj.% f 

s=. ~kC M-e“~^-e)C 

C “*^c^ + C W. ~ Ue)( U^. -i-He^ 

“ 0 + i T’iC^'> t 'P. + v"3 


^■pTe 


But 






d^Te 


T^ 


u|- ' 

where M is the jet Mach number. It folloxjs that 

; , s = i . 

° ,V-,-! ,1 - . c 6 - O + t. *■ '~^'} <e 


Solving for T^/T^j it is found that 

X 

Te 

A. 2 . Entrainment 

yf\ = 2Tr ^ surdr 


^ - 1 -^XOO = r & - 1 + lyU-D *■ 




= 2.Tr| 5’§urdr -5- 5^ 

^ surdr + is,iur’-- -L?,UeTf ] 


CA.35) 

The entrainment is the increase in the mass flow rate per unit distance 
Hence, 

o r’^‘ 

“ En'i'ratnyne'ni = 2-Tr^j CSu. S’^.Ue^Ti^^^Tf' 

VJt p 

T.et . 

(.1 


CA.36) 


J §U. -- Se O U5 f 5’ C S a ~ i /t) 'I. d 

•e 




(A. 37) 


^ "' I '* 

where R is the average gas constant for the misture and may be assumed 

to be R„ , and 

r - \Lli±j£2. i? 

-I 1+ ^ 

^ + 0 . 0 ^ 06 lC -t- o-op35-ni (A. 38) 

1 + 0.78'1'fS'P 
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The factor 


, 1 . 


Te'K 


in Eq- (A. 37) can be eliminated by using Eq. (A. 9). 


Hence, 


ii 

_ 2.Tr-^ ^ f Tt F, Fa ~\ 
>^T, F, 'VV,’' F. / 


(A. 39) 


(A. 40) 


According to Ricou and Spalding (Ref, 23) and Hill (Ref. 24), 
^r^Y":?2rL = consW -for 0 . 

Therefore, Eq. (A. 40) can be written as 


Ap r -;y>n „ — T,F 

*■ Wj. Uj' {%-f- 


But 


M 


_ _ J_ 


"^^o" 5jUj- U,- 

It follows that 


’•Tr 




^o( J J_ 2- 

\ sJ /«-T,F, + T,^Fx 

= 3:il2'-i} -5 T.F, 

4% ■)-?,’• F^ 

_ a. (l-y^> S t/F, •^ T, f/ -/^Tx Fa - 

7X ^ + T^’-Fx 


(A. 41) 


where ' — 




C?.n-/"T.F3)0^T^>,V*T^nV=^Ti'P/r=.+T^"^^^^ 

c/«T>.p,.T,*F.y ^ 

In consistence with the turbulent jet model, the 


(A.43) 


derivatives in Eq, (A. 42) will be evaluated far downstream and then 

eKtrapolated to the near fiead. Hence, 

j d r 1“/*^ ^ 

XilO< 5 ~ 1 -1- ic - 

- “ T, ® t "Fi 

■Pi C|^.) - 1 -h :L ^ 8} - 


T,' e F T, -b ^ 'Fi ® 'F.h) (A. 44 ) 


38 





L- -> ~ 2 H, 
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(A. 45) 


F,' = -T'V 411 li& 

” C 1 + T:.0-1^)} ^ 3 

at ^1 P-°?‘?o| - o.o4ooSTg+ o.ol (A.4 o) 

}+ l.os-ool?_ 

* ^ 

^ _ -p » 0.05-3.70 - — P. 0 ' 2 . -f *P^ ( A . 47 ) 

\ -^ Uog^ 6^ 

f^: 

' -i 1 ; ,+ T 

^ _.pt 0 - 0.0 4^A3Ti + o* °I S 2 . oT-..- ‘^•°'’ (A. 48) 

1 + I. 0Z2-72. Ti. 

One special application of Eq. (A. 42) is the isothermal mixing 
so that Pj = 0. Eq. (A. 42) is reduced to 


A 


’ A,S^22!!L — \ Ti* Fi (TiFt^C/^Ti fi -h^TtTi Fn.} 7 

•'• §e^ " 7 ^ lyU?,F,+?"F, C/^T,Fj-i-T7"F0" ^ 

/!T (./»?, 

C/^F,+ T,F.J . Vy- " • a 


= ^-^.FiF. 


T,' 


Ti 


c/^f.+t.f.T" 


This expression is independent of density ratio is k, is also inde- 
pendent of it. As ^ -5> o , the entrainment at the end of the potential 
core can be obtained: 


S-f 2IL = ( _L+^_{ 
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A special case is ^ = 0 


It can be shown from Eq. (A. 50) that the 


entrainment function is a constant and is equal to 


_£c 




■^x 


(A. 51) 


P. o 


which depends on the value of kj . Therefore, to complete the theory, 
it J.S necessary to choose k, such that best correlation with experi- 
mental data can be obtained. 

It is of interest to note that different values of k, have been 
used by different authors in the past (Refs. 25 and 26). In Prandtl’s 
model for uniform density mixing, it is assumed that 

£ = *, ( «.. - wo 

and k, has been taken as 0.022. In Ferri's model, 

§1 = C UO 

and k| was chosen as 0.025. In correlation with experimental data of 
free jets (with quiescent surrounding), Kleinstein found that k, = 0.0185. 
Recently, Witze (Ref. 27) re-correlated free-jet data and found that 

(A. 54) 

Unfortunately, if Eq. (A. 54) is used in Eq. (A. 51), it would result in 
the entrainment function depending on the density ratio, contrary to 
the experimental conclusion by Ricou and Spalding (Ref. 23). A more 
complicated model was proposed by Boyle and Viets (Ref. 25); 

(A. 55) 

where K = 0.022 in the main region of the jet. According to this model, 
the effective mixing length in the far field is greater for jets in a 
coflowing stream than for free jets because of the last factor in Eq. 

(A. 55). Comparing (A. 53) with Eq. (A. 55), it is reasonable to assume 
that k, in Eq. (A. 53) increases with the velocity ratio tAe/U| . 
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To keep the method as simple as possible, it is assumed in the present 
method that k, = 0.0185 for/* = 0, as used by laeinstain. At >* = 0.5, 
different values of kj have been tried and the results for the centerline 
velocity decay compared with experimental data. Reasonably good agree- 
ment lias obtained if k, = 0.024. This is shown in Fig. A. 2. For other 
values of /* , a linear interpolation or extrapolation gives 

•ft, = 0.0 1 S5- + 0.0 11/1 (A.56) 

To show the accuracy of Eq. (A.56), comparison with experimental data 
and the computation of Boyle and Viets is presented in Figs. A.l - A. 9. 

The turbulent Prandtl number is assumed constant and taken to be 0.72. 

It is seen that the results of the present method are reasonably good, 
except Fig. (A. 5a) where the experimental data show different trend 
from those of Fig. A. 2, even though the density ratio is only slightly 
different. 

With Eq. (A.56), Eq. (A.51) shows that for /a = 0.285. 

This is to be compared with the experimental value of 0.32 by Ricou and 
Spalding (Ref. 23) and Hill (Ref. 24). This difference may be attributed 
to the assumed form of 5K , as has been discussed by Ricou and Spalding. 

To see how the asstamed velocity profiles would affect the entrainment, 
let 

= !(; 1 + WSV5,) (A.57) 

Hc-Ue. 

Eq. (A.57) was used by Squire and Trouncer (Ref. 6). If Eq. (A.57) was 

used, Eq. (A.51) would become, for = 0, 6 = 1» - 1> 

(A. 58) 

%/ KX 

The ratio F| /F,. would be 1.725, instead of 1.926 in Eq. (A.51). This 
would further decrease the value of the entrainment function. Since 


Fj is proportional to (u - s while F^ is prc^jrtional to (u - u^) , 
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Figure A,1 Theory and Experiments for Turbulent Jet Mixing, 
p /p.-=l, T./T =1, !J=0.25 and M.=0 



©mo Experiments, talteu from Fig. 4 
of Boyle and Viets 

Boyle and Viets (Ref. 25) 

Prandtl Model (Boyle and Viets) 

— - — Present Method 


Figure A, 2 Theory and Experiments for Turbulent Jet Mixing 
P„/P.=l, T,/T =1, y=0.5 and M =0 
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Figure A. 4 Theory and Experiments 

P^/Pj=2, T^/T^=2, p- ( 
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for Turbulent Jet Mixing. 
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O Velocity ■? Experiments , taken from 
Temperature-J Fig. 15 of Boyle and Viets 

~ TempSSurel Viets (Ref . 25) 

Velocity -j Method 
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Velocity '! Experiments, taken from 
leraturei Fig. 13 of Boyle and Viets 
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Figure A. 5 Theory and Experiments for Turbulent Jet Mixing 

(a) p^/p^ = 1.19, TpT^=1.19, P=Q.5 andMpO.l 

(b) p /p,=1.296, T./T =1.296, p=0,25 and M =0.1 
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Figure A. 9 Theory and Experiments for Turbulent Jet Mixing 
p^/Pj=1.09, T^/T^=l., y=0.2 and M^=0.2 




it seems reasonable to assume that F, should be replaced by SF| and Fx 
by , where S may be called the shape factor of the velocity profile. 

If this assumption is made in Eq. (A. 51), then 

S3 222L g t-^Sl 

fj o,o6£.76S 

With k| = 0.0185, this value is to be equated to 0.32 to give S = 0.8907. 

The reasonableness of the shape factor can be demonstrated by using Squire 
and Trouncer’s velocity profile. If Eq. (A. 57) is used, F, = 0.14868. 

This can be reduced to the present value of 0. 12857 by using a shape 
factor of 0.865. Then F^ = 0.08618 by Eq. (A. 57) will be reduced to 
(0.865)^ X 0.08618 = 0.0645 which differs only 3% from the correct value 
of 0.06676. It may be summarized that for the present application, in 
all expressions developed so far, 

F| te ty o,8*?07 P, ^ a.'hd Fx t’y (o.S'fO?) Fa. . 

Fj remains unchanged; i.e. , the temperature profile id-ll be assumed to 

be correct. 

Entrainment in the initial region 

The entrainment function will be assumed to vary linearly in the 
initial region. At the end of the potential core region, the value of 
the entrainment function is obtained by the above theory. At the exit, 
it will be obtained by Abramovich’s theory. Since the latter does not 
account for compressibility effects, it will be applied in the following 
manner For a free jet xri-th uniform density and isothermal mixing- the 
entriaament at the jet exit can be obtained by Abramovich theory. Under 
other conditions, the entrainment at the exit will be changed proportion- 
ally in the same way as that at the end of the potential core being changed 
relative to the free jet value. 
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According to Abramovich’s theory, the inflow velocity at the exit 


for uniform density and isothermal mixing in a free jet is given by 
[p. 306 of Ref. 21] 

(A. 60) 


^ = 0.^7 C B-A) 


where A = 0, 316 and B = 0.450 (p. 282 of Ref. 21). Hence, 

=5 2. 7r $ V. 

4 r „(2-rp)2.'rrr;? fin 

^ - T m. -tt r 7~' ^ - c.. M 6) ^ o.\^^ (A.61) 

A. 3 Comparison with Knovm Results 

It should be noted that the present method is applicable to non- 


uniform density and/or non- i so thermal mixing with compressibility effects. 
The compressibility effects are accounted for through Kleinstein's theory. 

The model of Squire and Trouncer in uniform density and isothermal 
mixing is compared with the present one in Fig. A. 10. For a free jet, the 
Squire and Trouncer 's results are computed by using the data shox«i in 
Table 1 of Ref. 6; while the other two curves are constructed with the 
information given in Tables 2 and 3 of the same Reference. Note that 
to convert Squire and Trouncer ’s results to the present notation, the 


following relation is needed : 




(A. 62) 


Since hij =Trr<,*' > it follows that for - ?3 , 

r ^3 2H - (A. 63) 

It is seen from Fig. A. 10 that Squire and Trouncer ’s model gives 


higher entrainment in the initial region of the jet and lower values in 
the fully developed portion. The difference with the present method be- 


comes small as _/i is increased. 

To find the effects of entrainment on the aerodynamic characteristics, 
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Hill's Experiment (Ref. 24) 

Present Model 

Squire & Trouncer’s Model 



n of Entrainment Models, 
nd M.=0 


p^/p^=i. 



a line sink distribution is placed on the jet axis with strength equal 

to the entrained voJuae flow rate: 

o I (A. 64) 




Let 


Then 


’’’jA %} 


(A. 65) 


The corresponding induced velocity components due to the sink distribu- 
tion can be computed as (Ref. 15) 




(A. 67) 






(A. 58) 


■»x 4Tr j 


[tx-l j"+ys"r"f^ 


(A.69) 


In application, the jet axis is divided into segments within which |-C^) 
or EC|) is assumed to vary linearly. Let 

Eq)= 

in the ith interval. Then 

Zi- - _ J_ ^ U* r f ? (flj + W!)^! 

' _ _ C Sz)' ^ T" ^ C '^ !) 


or. 


(A. 70) 


Similarly, 

J_ Z± = r 

' d I ^ 

-’X. 




1 

i 3x. tcx-!)"+<a^f 




(A.71) 


Ths intsgirals in E<^s* (A. 70) and (A»71) can bs exactly integratsd to 
give 


^ =. - _Zii^ G,. - W Gxi 


(A. 72) 
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-b b. L +/c«- ^vnt (A. 73) 

+ /cx-A*0’' + y®^?' 

-;■ (A.74) 

^ , = ■ ■ .! - * (A.75) 


The above formulation is now applied to find the pressure distri- 
bution due to entrainment on an infinite plane wall with a cold jet ex- 
hausting normally from it. The wall boundary condition can be satis- 
fied by using an image sink distribution. The results of using this 
image system is to cancel the velocity component at x = 0 (i.e., 

on the If. all) and the velocity 24_ on the wall is doubled over the value 
when the wall is non-existing. Using the present model, the induced velo- 
city on the wall is given by 


J- 


Ujr. "»r 


where 


x=o 


ft =0.1^5' , O.OIJS , 

^ _ k - , X =Z0 


(A.77a) 


^j.= 




X V 


(A.77b) 


For the Squire and Trouncer’s model, the initial region is divided 
into 4 segments in accordance with the results shoim on Table 1 of Ref. 6. 
The entrainment data are then computed by the method indicated there. 

The results are as shown below: 
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= o-/ 3 ’- 5‘4 + 0.017*76 ^ 5; X <• 2,77 

— c. 17741 -fo.oU^S ^ , 2.71 :S 5c < 4,1-§ (A. 78) 

= 0. /iTI 33 4 0. 015 7.^7 it , 445 C , o '^ 


= 0 . 24 *^ 


, X > 6.o«^ 


I _ _ 

iAjr„ -^f 


^ = -zx 


^5 - <3,J + r + 7 =S^i£= 1 1 

X=0 , 


The pressure coefficient expressed in terms of jet dynamic pressure is 
given by 

C =_(J_ Ji-Y 

'®r ) (A. 80 ) 

The results are compared in Fig. A. 11. The experiment was done 
by Gentry and Margason (Ref. 28) and Wygnanski’s analytical results are 
given in Ref. 29. Both curves were taken directly from Fig. 1 of Ref. 
(30). It is seen that higher entrainment in the initial region in Squire 
and Trouncer's model yields slightly higher than the experimental re- 
sults and the present method. On the other hand, the present results are 
generally lower than the experimental values away from the jet boundary. 

It should be noted that small values of C are being compared which are 

P 

very sensitive to the change of entrainment values near the jet exit. 

Away from the wall, the inflow velocity can be computed as follows. 
Using the image sink distribution, it is obtained that 
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2^ 

'if 


^ { [cs-f 


^ zjl 5 r Ji 

47r 

oi:; ^‘^'> 


„_£l 

-fir 


where 


c'-i:: 

<= 


, ^0,5 2. , f>rese>ii- ■h\o<iel 

0.2.4‘J 
R-V Xc 


4*} I Ssjuire a~t.<4 TroiAVieey 's inode I . 
“ K 




(A. 82a) 




yiTT^ 


(A. 82b) 


The results on vertical planes normal to the jet axis and at x = 10 
and X == 40, respectively, are shown in Figs. A. 12 and A. 13. As expected, 
the lower entrainment in the far field in the Squire and Trouncer's 
model results in lower inflow velocity. 
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Appendix B 


Properties of an Equivalent Inviscid Jet 
A turbulent jet has its properties varying not only in the axial 
direction but also in the radial direction. The perturbation in such a 
nonhomogeneous flow region would be governed by a partial differential 
equation with non-constant coefficients. To simplify the matter, an 
equivalent inviscid jet with uniform properties is used in the inter- 
action computation in the present investigation. Such an equivalent jet 
is defined in accordance with the conservation of momentum, mass and heat 
content. Let the equivalent jet properties be denoted by ?. , Uj and*?- . 
It follows from these conservation principles that 


\ 2. Uo A r. / 


He 


(B.l) 






(B.2) 


O ^ 

TO- 


The left-hand side of Eq. (B.l) can be written as 


% Ue 

liHi ( 2lt-7r( ? 4 'F'F 

%v.i\ rj Z' V t ^ '^1 T^r 

.^here Eq. (A. 9) has been used . Hence, Eq. (B.l) becomes 

+ •n’-F-.)= i f|(-% - Oi\T (=•» 

K ) “ Similarly, the left-hand side of Eq. (B.2) 

can be written as 

U-|rj = 


2L^,1 r p p p . 
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where Eq. (A. 37) has been used. It follows that Eq. (B.2) becomes 

' /‘■R. - 0(|T 

To integrate Eq. (B.3), note that ^ 

i'u, Vr./ v S, -»!, 

= C-t)'|7r5 » \ 



(B.8) 

Hence, Eq. (B.3) becomes 


/if n-Tt = 


To solve Eqs. (B.5), (B.7) and (B.9) for the 

equivalent jet 

properties, let 


= +'R’'FO 

(B.IO) 


(B.ll) 

h = - ( t,f, -.^?.fo3 

(B.12) 

IV 

*,= JililL 

•Se Ue 

(B.13) 

Ug 

(B.l^) 

= c 

(B.15) 

Then, 



(71.16) 
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~ 0 

Solving Eqs. (B.16) - (B.18), it can be shown that 


(B.17) 

(B.18) 




fl 

5a.+ "53 




X X 

■»3 f 



(B.19) 


(B.20) 


(B.21) 

Xj — ' 

For non-heated jets, T. = T . In this case, Eq. (B-3) is not needed. 

1 ^ 

The above solutions, Eqs. (B.19) - (B.21), are still applicable if 
is equated to zero. 
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Appendix C 

Applications of Shanks’ Transformation to Evaluation of Series 


Consider the evaluation of the integral: 

F .. (c.i) 

’ i I -t TV. C ’ ‘ 

Assuming small , the following series can be obtained: 

•X. ^3 

= £>. US^T! 4 ~ o.oS^on 


+ ^ ° ^ i U ) 7 ^ 

By Shanks’ transformation, F, can be rewritten as 


F. S 


'hfl 


(C.2) 


(C.3) 


^ii+t Sii-i " ^ 


V 7 here is the sura of the first n terras in Eq. (C.2). Since (n+1) 
terms are needed, n will be taken as 4 in Eq. (C.2). Therefore, 

^hfl o, f 3. 7li- - o. oS o 1/ h '’■“667-S'?? Tj. o,oS-z-7cT-Ti +-o. o-f | T, 
After simplification, Eq. (C.3) gives 


o. i- p.ol6 )7 TV, + ^ T,- 

^ o. S i ? I 7 


Eq. (C.I) can be integrated exactly to give ^ 

^ i Tx *• 6 •P.T^ 1 + p/"+ T/^ V 


_JL_‘7? T> _ J_ (C,5) 

T/* J J > ' Tv 

Using Eq. (C.5) as the standard for comparis^. , Eq. (C.2) gives error 
of 22% while only 0.2% for Eq. (C.4) with = 1.0. 
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